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Abstract 

I present a simple numerical model based on iteratively updating sub¬ 
groups of a population, individually modeled by nonnegative real nnm- 
bers, by a constant decay factor; however, at each iteration, one group 
is selected to instead be updated by a constant growth factor. I discover 
a relationship between these variables and their respective probabilities 
for a given subgroup, summarized as the variable c. When c > 1, the 
subgroup is found to tend towards behaviors reminiscent of antifragility; 
when at least one subgroup of the population has c > 1, the population as 
a whole tends towards significantly higher probabilities of “living forever,” 
although it may first suffer a drop in population size as less robust, fragile 
subgroups “die off.” 

In concluding, I discuss the limitations and ethics of such a model, 
notably the implications of when an upper limit is placed on the growth 
constant, requiring a population to facilitate an increase in the decay 
factor to lessen the impact of periods of failure. 


1 Introduction 

Antifragilty is a growing area of research in complex systems, and classical 
examples such as the hydra-who is strengthen by stress, up to a point, simply 
by growing two new heads each time an old one is cut off-pervade the literature. 
My aim in this work is to provide as simple a numerical model as possible that 
exhibits the nonlinear overcompensations expected of an antifragile system. 

The concept is straightforward: a population, say of rabbits, is composed of 
several subgroups, say of different colors. At each time step, such as generations, 
the population as a whole is exposed to a stressor from a set of possible stressors, 
such as a potential predator. Each subgroup will decrease in number except for 
at most one-this subgroup may be well adapted to this stressor and actually 
benefit from its presence, increasing in number. 
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As time goes on, the population size as a whole will decrease until, if any 
such subgroups exist, the growth of the more often well adapted subgroups will 
offset the loss of the less adapted subgroups; eventually, it is hoped, the growth 
of these well adapted groups will overcompensate for the others, leading the 
population size to again grow as a whole. This is not unlike the self-healing 
material considered in [5], specifically those that “borrow from areas of less 
stress to fortify areas under more stress.” 

2 Model 

I define a population P as a set of positive real numbers H. By Hi I denote the 
ith subgroup of the population, and by the population size I refer to the sum of 
the elements oi H. 1 define a stressor 5” as a set of random variables such that 
St is a random group label i at discrete time step t. Belatedly, I define a set of 
random variables X such that Xu = /3 iff 5't = f; otherwise, Xu = a. By a and 
P I denote, respectively, a decay factor constant in (0,1) and a growth factor 
constant in (l,oo). 

The model proceeds as follows: at each time step, one subgroup is selected of 
the population at random (via St); the selected subgroup grows by the factor /3, 
and the other groups each decay by the factor a; if a subgroup ever goes below 
1.0, it is considered to have “died off” and it is set to 0.0; if the population 
size goes below 1.0, then all subgroups must have died off, and the model halts; 
otherwise, it continues until a predetermined number of time steps. 

Note that no assumptions are made about the probability distributions of 
St and, subsequently, Xu, other than that all St are assumed to follow the same 
distribution. This means that some subgroups may be selected more often than 
others. 

With this model, 1 am interested in the longterm behavior of the population 
size with respects to a given stressor: iff the population is fragile to that stressor, 
then it will tend to die off; iff robust, it will tend to an equilibrium; and iff 
antifragile, it will tend to grow infinitely, where an antifragile system is loosely 
defined as “a system that becomes stronger when stressed” [5], and in this 
particular application, a group exhibiting a nonlinear increase in group size 
over time in response to a stressor, i.e., set of random group labels. 


3 Analysis 

Let Hi{t) denote Hi at time step t, that is, after t iterations of the model; 
similarly, let P{t) denote the population size at time step t. Then, by (1)- 
(4), the expected value of Hi{t) is found to be an exponential growth or decay 
function, dependant on whether the geometric mean of Xi is, respectively, in 
(l,oo) or (0,1). The expected value of P{t), by (5), is simply a sum of such 
functions. 
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H, (t) = a”/3'"iLi (0), n -k m = t 

Definition 

(1) 

= (a"/*/3™/*)‘iL,(0) 


(2) 

lim H,{t) = (a^”[^*=“l/3^”[^‘=^l)*iLi(0) 

t—>-oo 

Law of Large Numbers 

(3) 

E[Hi{t)] = G[XifHi{0),G[-] = geo. mean 

Expected Value 

(4) 

E[P{t)] = J:,E[H,{t)] 


(5) 


Let /3 be chosen to be a function of a, the distribution of Xi, and a constant 
c such that (3 = ^nd c > 0. It is found by (6)-(9) that 

the fragility of a subgroup can be determined solely by its value for c in this 
function for (3. 


H^{t) 
c < 1 
c = 1 
c > 1 



= a"-^”iL,(0) 

Substitution from (2) 

(6) 

H,it) < H,{0) = 

G[X,] < 1 

Fragile 

(7) 


G[X,] = 1 

Robust 

(8) 

H,{t) > H,{Q) = 

G[X,] > 1 

Antifragile 

(9) 


This implies that E[Hi{t)] for each fragile subgroup t is a decay function 
and will tend towards zero, eventually dying off as it crosses below 1.0; and 
that E[Hi{t)] for each antifragile subgroup z is a growth function and will tend 
towards oo. By this approach though, E[E[i{t)] for each robust subgroup z is 
a straight line, implying that both robust and antifragile configurations will 
“live forever.” However, due to the model’s high variance, any subgroup can 
experience a long sequence of as, the population dying off as a whole out of 
“bad luck.” 

So, let L[-] denote the lifespan of a variable, equivalently, the expected num¬ 
ber of time steps before the variable first falls below 1.0. Next, for a given 
subgroup z, consider the sequence of as and j3s in X^, ordered by t. If we define 
Xit for time steps only during which z has not yet died off, then obviously not 
all possible orderings of as and /3s are possible; for example, a sequence of n as 
followed by a single (3 is not possible when a”iLi(0) < 1.0 because the subgroup 
would have already died off before reaching the /?. 

Let K be the random variable of the number of /3s in X^ for a given subgroup 
i. Let the initial health constant w = logQ,(iL~^)-|-l represent the initial number 
of as needed for subgroup z to die off and let the compensatory health constant 
w' = logQ,(/3“^) represent the number of extra as needed for that subgroup to 
die off after that subgroup “sees” one /3. 

Therefore, L[E[i\ = w -\- E[K]w' -\- E[K], that is, the sum of the initial health, 
compensatory health for each /3 seen, and the number of /3s seen themselves. 
By (10)-(15), the probability mass function of K is found in a manner similar 
to that of the binomial distribution. 


3 



L[Hi] e {w + kw'\k = 0,1,...} (10) 

Let K = {L[Hi] — w)/w' = k (11) 

Pr[K = 0] = Pr[Xi = a]™ (12) 

Pr[K = k] = Pr[Xi = a] = l3fC{k) (13) 

Where C{k) = valid orderings of w + kw' as and k /3s (14) 


Figure 1 examines the commulative distribution function (CDF) of K under 
several different configurations, illustrating that the behavior of the distribution 
is heavily influenced by c: when c > 1, that is the subgroup is antifragile, the 
CDF quickly converges to a value less than 1.0; and when c < 1, that is the 
subgroup is fragile, the CDF quickly converges to 1.0. This general behavior is 
irrespective of Pr[Xi = a] and w. 

Because Pr[K = k] represents the probability that subgroup i will see ex¬ 
actly k /3s and then die off and Pr[K < oo] represents the probability that the 
subgroup will die off at all, the behavior of the CDF implies that an antifragile 
subgroup has a significant probability of “living forever,” whereas fragile sub¬ 
groups have almost none and robust subgroups are not as straightforward to 
predict. 


4 Verification 

To verify the predictive power of the model, 5,000 simulations were run for each 
k from 0 to 100 where /3 = 2, subgroups were labeled 1..5, subgroup i was se¬ 
lected for growth ^ times, and iLi(O) = 10 for each subgroup. The remaining 
parameters are found in (16) - (24) such that subgroup 3 is expected to be 
robust, subgroups 1 and 2 are expected to be fragile, and subgroups 4 and 5 
are expected to be antifragile. Figure 2 illustrates the results, comparing cal¬ 
culated (predicted) CDFs for each subgroup and observed (simulated) survival 
rates. Predictions were made in a bignumber implementation of Julia 0.3.5 and 
simulations were run on a 64-bit floating point implementation in Python 2.7.8. 
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Antifragile(c=1.75) 


Fragile (c=0.25) 



-Pr[a]=0,1 

.Pr[a]=0,2 

-Pr[a]=0.3 

— Pr[a]=0.4 

-Pr[a]=0.5 

-Pr[a]=0.6 

-Pr[a]=0.7 

.Pr[a]=0.8 

-Pr[a]=0,9 


Robust (c=1.0) 



Figure 1: CDFs of several Ks: antifragile (top left), fragile (top right), and 
robust (bottom), with Pr[Xi = a] from 10% to 90%, k from 0 to 100, w = 5, 
and w' = 1 for each. 
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p 

= Q,-'P'-[Y=a]/Pr[Js:i=/3] 

Robust 

(16) 

p 

= 2,Pr[Xi = p]= i/15,i = 3 

Given 

(17) 

2 

= a~‘^ 


(18) 


^ a = 2"^/^ « 0.84089 


(19) 

w 

= logM(O)-^) + 1 

Definition 

(20) 

H,{0) 

= 10 

Given 

(21) 

w 

= 1032-1/4(10-^) -b 1 « 14.2877123 


(22) 

w' 

= logaiP~^) 

Definition 

(23) 

w' 

= ^0ff2-i/4(2-^) = 4 


(24) 


Figure 2 also illustrates the behavior of a population size over time for a 
model configuration with 100 subgroups, Ffi(O) = 10 for each subgroup, /? = 1.5, 
a = 0.995, and t < 200. I compared two methods for implementing St', an iter¬ 
ative method where St = {t mod 51) -I- (mod 51); and a random method 
where St = uniform(0, 50) -I- uniform(0,49). Note the clear global convexity of 
the population size over time-a requirement of antifragility-, although this con¬ 
vexity is composed of several local periods of concavity in the iterative example. 
In the random examples similar patterns can be noted, although they are not 
as pronounced. 

5 Conclusion 

I have presented and analyzed a simple model that exhibits antifragility. In this 
model, the population as a whole is exposed to a stressor, represented by a set 
of random variables. At each time step, exactly one subgroup of the population 
benehts from this stressor, growing in size as the others shrink. In the long 
run behavior of the model, the weaker (fragile) subgroups of the population 
die off, represented as an exponential decay function while the stronger (robust 
or antifragile) subgroups either maintain an approximate equilibrium or grow 
infinitely in size. 

The model either contains no or masks with its simplicity any feedback 
loops, posited in [1] as necessary to produce a stable system: “positive feedback 
alone pushes the system beyond its limits and, eventually out of control, while 
negative feedback alone prevents the system from reaching its optimal behavior,” 
a behavior echoed here in the tendancy of fragile groups to die off, antifragile 
groups to grow “infinitely,” and robust groups to rest somewhere in between. 
Without loops of any kind, then, the proposed model is instead little more than 
a set of functions with tendancies to grow or decay in response to a set of random 
variables. 

Yet does this model still suggest brute survival of the fittest? Yes, but only 
in cases where it is directly applicable, and these applications presuppose a 
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Predictions vs. Results 



k 


Iterative 


Stochastic 




Figure 2: Comparison (top) of predicted CDFs (lines) and observed survival 
rates (points) for subgroups 1 (top line) through 5 (bottom line). Also, popula¬ 
tion sizes over time using an iterative method to select subgroups (bottom left) 
and a random method (bottom right). The iterative data is shown on the right 
alongside three random runs’ data for quick comparison. 
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static longterm stress environment and complete independence of population 
subgroups, neither of which are typical givens in complex systems such as so¬ 
cial networks and biological ecosystems. This is, as it has been characterized 
here, the model does not support a stressor that changes suddenly, presenting 
an opportunity for the model to suffer the “turkey fallacy” [3] , where the system 
is fragile to abrupt environmental changes. Furthermore, because of the afore¬ 
mentioned limitations on the model, a general measure of antifragility cannot 
be given, albeit a general sense of the term is. 

However, it does demonstrate a relationship between periods of failure, repre¬ 
sented by the decay constant a, and periods of overcompensation, representated 
by the growth constant /?. If an upper limit is placed on the ability of a popu¬ 
lation to overcompensate from a stressor, as one would expect is often the case 
in natural systems, then once that limit has been reached the population’s only 
choice to improve its robustness is to raise the decay factor-that is, increase a 
through some facility such that during periods of failure when overcompensation 
is not possible the impact of the stressor is not as severe. 
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